Renormalization Group and Curved Spacetime 



H. Matsueda 

Sendai National College of Technology, Sendai 989-3128, 
(Dated: June 29, 2011) 



Japan 



o 

(N 
C 

oo 



o 



in 



i 

C 

O 

o 



> 

(N 

in 

O 



X 



We examine the role of curved geometry on renormalization group by means of image compression 
based on the singular value decomposition. By calculating course-grained images and their entan- 
glement entropy, we find the anti-de Sitter space / conformal field theory correspondence hidden 
in the compression. The correspondence is originated from the conservation law for information. 
We discuss how one particular metric is automatically induced from various images. A formula for 
resolution of the course-grained images is also derived. 
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The renormalization group (RG) or the course graining 
is a fundamental concept in theoretical physics. Renor- 
malization was originally considered to be a method for 
avoiding difficulties of divergences in quantum field the- 
ory, and RG was applied to critical phenomena in con- 
densed matter physics. However, we may not still get to 
its heart, since the recent progress of RG requires more 
comprehensive understanding among entanglement, du- 
ality, holography, curved spacetime, and so on. 

A key concept concerned with renormalization is the 
anti-de Sitter space / conformal field theory (AdS/CFT) 
correspondence [ij-yj]- This is holographic correspon- 
dence between a quantum field theory in spatially d- 
dimension, CFT^+i, and general relativity on hyperbolic 
geometry in one-higher dimension, AdS^+2- Thus, the 
correspondence is closely related to quantum gravity, and 
at the same time it has been shown that the radial axis 
of AdS space corresponds to a renormalization flow pa- 
rameter for CFTd + i that is located at the boundary of 
noncompact AdS^+2 manifold The role of this RG 
flow on course-grained information is reflected in the en- 
tanglement entropy. The entropy of a subsystem in a 
CFT is proportional to the area of the minimal surface 
in AdS space surrounding the subsystem [f| . This means 
that the total amount of the original information in the 
CFT is distributed to each RG layer. 

In condensed matter physics, the real-space RG is a 
long-term topic. There are two main streams: one is 
Wilson's pioneering work for the Kondo problem Q, 
and the other is White's density matrix renormaliza- 
tion group method (DMRG) 0- Recently, it has been 
found that the Wilson's RG controls the energy scale of 
low-lying excitations by exponentially bending the model 
space [1] . A systematic method for this energy-scale con- 
trol is called hyperbolic deformation [9]. On the other 
hand, DMRG is known to be a variational optimization 
method of the matrix product state (MPS) for gapped 
systems, and its extention to higher-dimensional systems 
is given by the tensor product state or the projected en- 
tangled pair state. The application of tensor networks to 
critical systems leads to the tree tensor network (TTN) 
and the multiscale entanglement renormalization anzats 
(MERA) [H[ll|. They are spatially (d + l)-dimensional 
hierarchical networks that are compatible with real-space 



RG, and the network forms discrete AdS space. A crucial 
improvement in MERA is to introduce the disentangler 
transformation. In a viewpoint of the AdS/CFT corre- 
spondence, however, the network structure seems to play 
a role on succesive description of critical phenomena [12j | . 
Therefore, the hyperbolic geometry is a common ingre- 
dient in the two different RGs. 

In order to explore the role of curved geometry on RG, 
we examine image compression based on the singular 
value decomposition (SVD) that is a core algorithm of 
DMRG and principle component analysis. Our approach 
is advantageous to extracting geometric information from 
images, since the color degrees of freedom can be re- 
garded as the strength of the gravitational field. Fur- 
thermore, various field patterns are generated by only 
exchanging the images. We find the AdS/CFT corre- 
spondence hidden in the compression of any images. The 
entanglement entropy of course-grained images shows the 
universal scaling behavior in one-dimensional (ID) quan- 
tum systems at criticality. On the other hand, the course- 
grained image itself is devided into patches represented 
by linearly independent bases. Then, the RG flow of the 
patches forms discrete AdS space. We have a conjecture 
that the independency is related to the conformal sym- 
metry. Based on the holographic entanglement entropy, 
we derive a formula for resolution of the course-grained 
images. The formula suggests the origin of the corre- 
spondence that comes from the conservation law for the 
amount of information. 

Let us start with performing SVD to image data 
with the pixel size M x N. We restrict our attention 
to grayscale images, but its extention to RGB data is 
straightforward. Then, the image data at pixel {x,y), 
^(x,y), is represented by an integer ranging from to 
255. A white (black) pixel is assigned to be 255 (0). We 
regard ^>(x,y) as a matrix, and ^S(x,y) is normalizaed 
to be ^(x,y) oc ip(x,y) so that J2 x , y ^ 2 ( x >y) = 1 - % 
applying the SVD to ip(x,y), we obtain 



1>{x,v) =Y t Ui(x) y /\ l Vi(y), 



(1) 
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with use of sigular values, A;, column unitary matrices, 
Ui(x) and Vi(y), and L — min(M, V). We assume Ai > 



2 



A2 > • • • > Al without loss of generality. 

The singular values are also the eigenvalues of the re- 
duced density matrices defined by 

P x(x,x') = '£il)(x,y)1>(a/,y) = ^U l (x)X l U l (a/){2) 
y i 

PY (y,y') = £>(z,2/)V>(z,2/) = X>(2/) A * W)-(3) 

v i 

In a viewpoint of DMRG, X- and ^-components of 
^){x,y) behave as system and environment blocks of ID 
systems. We thus expect that A; decays exponentially as 
I increases. Then, taking m-largest singular values is a 
good approximation for tp(x, y) (m <C L) 
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A quality of the course-grained image, ip m (x,y), is eval- 
uated by the entanglement entropy defined by 



S m = Xi logAj. 



(5) 



Note that S m picks up long range correlation more pre- 
cisely than the truncation number in DMRG defined by 
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FIG. 2: Zebra image, (a) m — 1, The black vertical line 
represents his nose and front leg. The lower half of this image 
shows oscillation of colors which comes from stripes of his 
body. The algorithm is going to replace the original image 
by one typical basis, (b) m = 2, (c) m = 3, (d) m = 4, (e) 
m — 20, (f) original image before compression (m = 425). 




A,,=0.859 (zebra) (c) 
A,,=0.874 (car) 
^.,=0.777 (flower) 
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FIG. 1: Sm for various images: (a) zebra (o, 425 x 640), 
car (A, 512 x 512), my office (v, 512 x 384), colleagues 
(•,512x384), cat (□, 398x342); (b) flower (x, 640x425), two 
cats (+,398 x 342), fractal tree (0, 1024 x 768), random dots 
(dashed line, 500 x 500), (c) X m for zebra, car, and flower. 

Figure Q] shows S m (and A m ) for various images as a 
function of log rn. Figures 2 and 3 present ?p m (x,y) for 
two different types of images (zebra and flower). Here, 
the pixel sizes of zebra and flower images are 425 x 640 



(horizontal x vertical) and 640 x 425, respectively. Then 
we have 425 independent eigenstates for both cases (L = 
425). In Fig. [T] (a), we clearly see that the various image 
data fit one straightline for small m region: 



1 , 

6 



(6) 



with 7 = — AilogAi ~ 0.13. Here, we assume that 
the images are in the open boundary condition. This 
scaling behavior is also seen in the MPS formulation 
(rn represents the matrix dimension in this case) of the 
transverse-field Ising chain at criticality (central charge 
c = 1/2), H = -Ei(°?°?+i + of) P> and the S = 
1 XX Z chain with the uniaxial anisotropy (c = 1), 
H = EiiSfSf+. + SfSV^ + ASfSt^ + Dj:^!) 2 , with 
A = 2.59 and D = 2.3 QJ. They belong to the different 
universality classes, and thus the eq. © is recognized to 
be universal for ID quantum critical systems. When we 
take larger m values, S m starts to deviate from the scal- 
ing line and approaches a saturated value. It depends on 
images how abrupt the saturation occurs. The saturated 
value tends to increase with M and AT, since the entropy 
measures the amount of information. Figure [2] (e) shows 
that the original image is almost recovered for m ~ 20 
(log m ~ 3) where the deviation starts. 

Figure Q] (b) shows exceptional cases that do not obey 
the eq. ©■ For the flower image (cross symbols), the 
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FIG. 3: Flower image, (a) m = 1, (b) m = 2, (c) m = 3, 
(d) m — 4, (e) m = 20, (f) original image before compression 
(m = 425). 



entropy for small m region has a gradient two times larger 
than the previous examples: 



S m = ^logm + 7, 



(7) 



where the residual entropy Si = 7 is still unchanged. 
According to the CFT 1+1 , the difference of the gradient 
in the entropy comes from boundary conditions. In the 
present case, the scale invariance among flowers seems 
to induce the change in the gradient (see Fig. [3]). For 
the tree image in Fig. [T] (b), the gradient of the entropy 
jumps from 1/6 to 1/3 at around m ~ 20. This comes 
from the fractal nature that the branches and the trunk 
have quite similar structures. The branch structure can 
not be seen for m < 20, while their fine structures ap- 
pear almost abruptly as we increase m. According to the 
Zamolodchikov's c-theorem, the fractal case is special, 
since the c-function corresponding to the gradient of S m 
does not change monotonically |15| . 

Let us next look at the course-grained images that pro- 
duce the eq. fl5J). We find that the image space is dev- 
ided into smaller patches as m increases. For instance 
in Fig. [3] (a) with m = 1 the image is characterized by 
one horizontal line, while in Fig. [3] (b) with m = 2 we 
clearly see discontinuity between upper and lower halves 
marked by a left arrow. The discontinuity comes from 
a fact that the upper and lower halves are represented 
by linearly independent basis states. This point may be 
formulated more precisely by the multi-resolution anal- 
ysis, and then we conjecture that the scale factor of the 
mother wavelet plays a role on the conformal symmetry. 
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FIG. 4: (a) Discrete AdS space, (b) Geodesic curve on conti- 
nous AdS space (CFT is defined on the a-axis). 



In the present case, the poistion of the partitioning is 
automatically selected so that the linear independency is 
optimized by DMRG. If we regard these patches as inde- 
pendent lattice sites, the lattice plus the RG axis forms 
the discrete AdS space shown schematically in Fig. [4] (a) . 
The construction of our AdS space is analogous to that 
of the MERA network with disentangler. TTN is also 
similar to the AdS space, but the tensors are still entan- 
gled. Therefore, the disentangling processes are crucial 
for generating the AdS space. 

Note that Fig. 2] (a) is not a oversimplifyed view to an- 
alyze the scaling relation even quantitatively. The conti- 
nous AdS metric has various expressions. For instance, 



ds 2 



_2 



(dz 2 + dx 2 ) = (dr log a) 2 + (a- T/l dx) 2 (8) 



with curvature I and z = ol t I 1 1. The factor a deter- 
mines how many patches merge together, but this trans- 
formation guarantees that the most simple expression, 
ds 2 = (l 2 /z 2 )(dz 2 + dx 2 ), is essential. 

Let us derive the resolution of the course-grained im- 
ages by calculating the holographic entanglement entropy 
within a continuous approximation (see Fig. |U(b)) @: 



S = ~~~ 



1_ 

4G" 



(9) 



where 7 represents the geodesic distance between X and 
Y on the same RG layer with m, and G is the Newton 
constant. The geodesic equation is given by (d 2 x x /dt 2 ) + 
T^(dx^/dt)(dx"/dt) = Owith x» = (x\x 2 ) = (x,z) and 
the Christoffcl symbol A solution is a half cycle: 

(x,z) — (L/2)(cos8,sm8), where the angle 9 obeys the 
following relation dO/dt = sin 8 (e < 9 < ir — e). We 
introduce a parameter a = (L/2) sine instead of e. The 
geodesic distance 7 is then calculated by 




I log 



1 + cos e 
1 — cos e 



(10) 



and the entropy is given by 



S = Tclog 
6 



L + v/i 2 - (2a) 2 )' 



(11) 
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FIG. 5: (a) A m and (b) S m for the zebra images with different 
pixel sizes: 425 x 640 pixels (open circles), 212 x 320 pixels 
(open triangles), and 85 x 128 pixels (crosses). 



where c = 3Z/2G [Tfj. Usually, we take the limit a — > 0+ 
in order to calculate the entropy of the original quantum 
systems. Then, we obtain S — (c/3) \og(L/a) [l7(- On 
the other hand, here we are going to approach 2a — > L, 
and the entropy should be comparable with the scaling 
relation ©. Thus, we obtain 



L+^/L 2 - (2a) 2 



> 2c 



2a 



(12) 



Finally, the resolution of a course-grained image for a 
fixed m value, R m , is given by 



Rn 




(13) 



This formula is consistent with numerical data, when we 
take c = 1. Here R 2 = 0.333, R w = 0.818, R 20 = 
0.905, and Rao = 0.951. According to this trend, we have 
already mentioned that the data with to = 20 almost 



reproduce the original one except for very fine details. We 
have confirmed that the compressed data with to = 10 
starts to reproduce the original one, but fine structures 
are still lost. The reasonable correpondence between the 
holographic theory and the numerical data also supports 
the existence of the AdS/CFT duality. 

Based on the duality of coupling parameters between 
AdS and CFT, we consider the reason why the Gaussian 
universality class, c = 1, appears. The zebra image repre- 
sents an 'ordered' pattern induced by strong interactions 
in this image space. Then, we can interplete the present 
results as a consequence of free bosons in the CFT side. 
We can exclude the possibility of taking another c values 
such as the Ising universality class c = 1/2 and other uni- 
tary discrete series, since the images do not have any Z 2 
symmetry. Furthermore, we always obtain the positive 
residual entropy 7, and this is due to weak topological 
nature of ID systems [ill, [l9j . When we look at the ze- 
bra image, he has various stripe patterns: The stripes of 
his body are wide, his face has somehow narrower stripes 
than his body, and his hair is more complicated. In par- 
ticular, his hair represents high frequencies, which corre- 
spond to low-energy behavior in the CFT. Actually, when 
we increase resolution of the zebra image by taking larger 
M and N values, fitting to the scaling relation becomes 
much better for larger to values (see Fig. [5]). This data 
clearly show the infrared cut-off in the CFT. 

Summarizing, the core of the AdS/CFT correspon- 
dence lies in the conservation law for information between 
the quantum and classical sides: The high probability 
bases of the reduced density matrix in the quantum side 
represent course-grained geometry in the classical side 
(see also Fig. [T](c)). The curvature of the classical space 
determines the amount of information that should remain 
after RG transformation. 
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